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The effect of interchain hopping on localization has been investigated for two coupled chains of 
spinless fermions with nearest neighbor interaction. By use of the renormalization group method, 
it is found that interchain hopping plays a significant role when the correlation gap induced by 
interplay of interchain hopping and intrachain interaction becomes larger than the Anderson 
localization energy scale induced by impurities. The strong competition between localization 
and derealization is examined for attractive interaction. 
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§1. Introduction 

The interplay of disorder and interactions in low di- 
mensionality is one of the most important topics in con- 
densed matter physics. In the recent years, it has known 
a renewed interest due to experimental evidences for a 
metal insulator transition in a two dimensional electron 
gas 1 ) in the absence of a magnetic field. Such a transi- 
tion being ruled out for non-interacting electrons, 2 ' the 
interplay of disorder and interaction is believed to be 
responsible for the transition. However, the theoreti- 
cal study of the interplay of disorder and interaction in 
two dimensions is an extremely difficult topic 3 ' and this 
problem remains extremely controversial. There exists 
also quasi-one dimensional systems in which the inter- 
play of disorder and interactions can be studied experi- 
mentally, such as quantum wires 4 ' and nanotubes. 5 ' In 
one dimension, the theoretical situation is better since 
bosonization techniques 6 ' can be used to deal with in- 
teractions, allowing for a Renormalization Group (RG) 
treatment of disorder. 7 ' It has been found that for a 
one-chain system, attractive interactions of the order of 
the bandwidth were needed to suppress Anderson local- 
ization and produce a metallic state. 7 ' More recently, 
it was found that in a model of two coupled spinless 
fermions chain 8,9 ' the opening of gaps in some excita- 
tions of the system could lead to a drastic modification 
of the response of the system to disorder. Namely, in the 
case when disorder is too small to close the gaps of the 
spinless fermion ladder, the system is delocalized for at- 
tractive interactions (see figure 1). As a result, for small 
enough disorder, the metal insulator transition occurs 
for weakly attractive interactions. However, the approx- 
imation of disorder weak enough not close the gaps of 
the system used in Refs. 8, 9 is not valid close to the 
non-interacting point where the gaps vanish, precluding 
a quantitative description of the metal insulator tran- 
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Fig. 1. The phase diagram derived in Ref. 9 For repulsive interac- 
tions, a localized insulating phase is obtained. For attractive in- 
teractions, a metallic phase is obtained. In the region V £ [V c , 0] 
the localization length is of the order of magnitude of the cor- 
relation length of the gapped excitations of the system. This 
transition region cannot be described with the methods 9 ' used 
to describe the two other regions 



sition in this regime within this approximation. Fortu- 
nately, the fact that the transition occurs for weak dis- 
order and weak interactions implies that the techniques 
of perturbative renormalization group can be applied to 
describe the transition. In the present paper, we under- 
take such an RG analysis. We will show that there exist 
three regimes. For sufficiently repulsive interactions and 
sufficiently large interchain hopping, there is a regime of 
pinned CDW n in which interchain hopping is relevant. 9 ' 
For attraction strong enough, there is a metallic phase in 
which t± is relevant. 9 ' Finally, for weak interactions or 
small interchain hopping, there is a regime in which the 
system behaves as a single chain insulator and interchain 
hopping becomes an irrelevant perturbation. 

The plan of the paper is as follows: In section 2 we 
recall the definition of the model. Then, we set up the 
bosonized Hamiltonian and the renormalization group 
equations. In section 3 we discuss the physics obtained 
from solving the RG equations. We identify the fixed 
points of the RG equations and discuss the nature of the 
corresponding phases of the two-chain system. Then, 
we determine the phase diagram resulting from the RG 
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equations. Finally, we discuss the behavior of the local- 
ization length with disorder strength. In section 4, we 
summarize our results and discuss possible extensions. 
Technical details have been left to Appendixes A and 
B. 

§2. Model and Renormalization Group Equa- 
tions 

2. 1 Hamiltonian 

We consider a system of two coupled chains of spinlcss 
fcrmions, described by the Hamiltonian: 
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where Cj j7 j (n=l,2) denotes a fermion operator in the n-th 
chain. The lattice constant is taken as unity. Quantities 
t and £j_ denote intrachain hopping energy and inter- 
chain hopping energy respectively ( t > t±) and V is 
the interaction between nearest neighbor sites. Quanti- 
ties £^1 and £^2 are the random potential acting on the 
chain 1 and the chain 2 respectively. One has: £,i lP £,i'y — 
W8i t ii6 Pi pi. In the following, we assume that t± is small 
enough that there arc two bands at the Fermi level. If p is 
the average number of fermions per site, this means that 
t± <C 2tsin (irp). We assume incommensurate filling. 

By diagonalizing the tj_-term and making use of the 
bosonization, eq.(2.1) is rewritten as 9,10 ' 
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In Eq.(2.2), random potentials expressed in terms of 
the symmetric and antisymmetric parts arc treated as 
UxjtaW) = D a S(x-x<) and|^X*(^) = D s 6(x-x<). 
We have neglected forward scattering since it does not 
contribute to Anderson localization. 11 ' The quantity 
m = — measures the difference of Fermi wavevec- 

u 

tors between the bonding and the antibonding band. 



The parameters K p , K\\ 
to t,tj_,V,W by: 
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where a is a lattice spacing and vp — 2ta sin(fcpa) is the 
Fermi velocity of a single chain. 9 ' The bosonized Hamil- 
tonian in the absence of disorder has been previously 
derived by Nersesyan et al. 10 ' by mapping the problem 
onto a Luttinger liquid with interactions breaking spin 
rotation symmetry in a magnetic field. 12 ' 

2.2 Renormalization Group Equations 

It will be convenient to write the RG equations in 
terms of the following dimcnsionlcss quantities: 
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In the following, we will take u p = uu = vp. By taking 
into account the effect of interchain hopping, the RG 
equations are derived as: (See Appendix A) 
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which are calculated with the initial conditions, 2//(0) = 
ia(0) = K p - 1 = -^sin(fc F a). For V a = V s = 0, 
these RG equations are identical to the ones of the pure 
system 
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2.3 The case i± = 

In the limit of t± = 0, a simplified form of the RG 
equations (2.8) can be obtained. With initial conditions 
#H (0) = 1, y/(0) = y ± (0), D o (0) = V s {0) and £j_ = 0, it 
is easily seen that one has K\\(l) = Ay(0) = 1, y/(l) — 
y±(l) and T) s (l) = D a (l), t±(l) = for any I. By using 
Vf — V± — K p — 1, it is seen that the RG equations (2.8) 
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can be reduced to a pair of RG equations for yj and 
T> = T> a ,s- the single chain RG equations are written as 

dl 

^ = {l-2y f )V . (2.10) 

These RG equations are identical to the RG equations 
for the single spinlcss fcrmions chain: 7 ' 

<IK„ .,..,.„ (2n) 



dl 
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in the weak interaction limit where K p ~ 1 + yf- As 
a result, for t± = we recover the RG equations of a 
single spinless fermion chain in the vicinity of the non- 
interacting point. Let us note that the scaling approach 
will not reproduce the prefactor in the localization length 
that is obtained in the case of non- interacting electrons. 14 ' 
From the limit t± = 0, we expect that for i± small 
enough, the two chain system will behave identically to 
a single chain system. For larger t±, a crossover to the 
ladder regime 9 ' should occur. 



with the metallic phase of the two chain system. 9 ' In 
the following subsections, we discuss in more details the 
different fixed points. In section 3.2, we discuss qualita- 
tively the transitions between these three phases. 
3.1.1 V{1) -► oo, y f (l) -► -oo, #|| (I) -► +oo 

This regime is represented on figure 2. Wc have K\\ — > 




2.4 The non-interacting limit 

As in the case of the single chain, the non-interacting 
limit has to be analyzed very carefully. It is known that 
in the single chain case, disorder seems to induce an effec- 
tive interaction in the system 7 ' even in the absence of in- 
teraction in the pure system. This artefact can be cured 
by separating properly disorder and inelastic processes 7 ' 
, so that starting from a non-interacting problem, one 
remains with a non-interacting problem. The effective 
interaction can be evaluated in the case of the two chain 



problem. If wc define as K 
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With these initial values, it can be shown easily that 
if we start from T> a = T> 8 and a non-interacting system, 
we will preserve Ku — 1, K p — 1, yf — y± = under the 
RG flow and thus no spurious interaction is generated. 

§3. Competition of disorder and interchain hop- 
ping 

3.1 Analysis of the RG flow 

In the present section, we discuss the strong coupling 
fixed points of the renormalization group equations (2.8). 
We show that there are three types of strong coupling 
fixed points. The first one is associated with the pinned 
CDW V of the two-chain system, 9 ' the second one is as- 
sociated with Anderson localization in the single chain 7 ' 
i. e. to an irrelevant t± and the last one is associated 



Fig. 2. The ^-dependence of t±(l), 
V s (l) and V a (l) for t ± (= t±{0)) = 
1CT 6 and K p {= K p (l)) = 0.9. 



K p (l), Kj|(J), y ± (l), y f (l), 
0.1,V(=V S (0)) = V a (0) = 



+00 so that flu becomes classical. Since j// — > —00, the 
expectation value of 6\\ is (6\\) = 0. This corresponds to 
the formation of the CDW n in the pure system, thus this 
regime corresponds to the pinned CDW n obtained for 
repulsive interactions. 9 ' Wc define the correlation length 
associated with this regime by: 1/£|| = c~'n (correlation 
length of 6*|| ) and l/£i oc . = c~ limp (localization length) 
where yf(l\\) = 1 and V a (l imp ) = 1. The energy scales 
corresponding to these correlation lengths are given by 
A|| = VfCT 1 ^ and E\ oc . = vfe~ limp , respectively. Within 
this regime, there are two possibilities. 

The first one, shown on Fig. 2 is that yf and T> a 
reach strong coupling simultaneously. This regime can be 
identified with the single chain regime by noting that Eq. 
(2.9) that describes the regime t± = has the same type 
of strong coupling fixed point. Let us note that a regime 
where T> a reaches strong coupling before yf is excluded 
since in Eqs. (2.8) T> a contributes to the renormalization 
of y f - 

The second possibility, shown on Fig. 3, corresponds 
to yf reaching strong coupling before T> a . In this regime, 
to obtain the localization length, one needs to renor- 
malizc up to the scale I* where yf(l*) = 1, then use 
the simplified renormalization group equations 9 ' that are 
valid when the expectation value of On is well defined. 
This regime corresponds obviously to the pinned CDW 7 ' 



regime. 9 ' 

3.1.2 2/ / (0->+oo,.K'||(0-»+oo,:D(0-> 
This regime is represented on Fig. 4. 
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Fig. 4. The ^-dependence of t±(l), K p (l), K\\(l), y±(l), Vf(l), 
Ds(l) and D a (l) for t± = 0.1, D = 1CT 7 and K p = 1.1. 



length ^i c h., the two chain localization length ^^° ano - 
the correlation length in the antisymmetric modes £y . 
It is clear that the two-chain localization length in the 
repulsive regime £2^° 1S always shorter than the single 
chain length £i c h.- Let us first consider the repulsive 
regime. For V, t± large enough and D weak enough, we 
expect £|| to be the shortest length in the problem. In 
this regime, the calculations of Ref. 8, 9 are valid and 



the localization length in the system is £ <~ £. 



v>o 

2ch. • 



This 



is the regime represented on Fig. 3. Let us now imag- 
ine that we are decreasing t± or V so as to increase £11 . 
We will first reach a regime where £11 ~ £2ch°- ^ n this 
regime, the calculations of Ref. 8, 9 are not valid any- 
more and the localization length £ starts to cross over to 
the single chain regime, i.e. ^X° < £ < £ lch -- This case 
corresponds to figure 2. Finally, as £11 ~ £i c h., the corre- 
lation length falls into the single chain regime £ = £i c h.- 
This is in particular the regime that obtains for t±_ = 
or V = 0. Suppose now that we start from V = and 
decrease V. Then, we are in the attractive regime, and 
£|| starts to decrease. When £11 ~ £i c h., we fall into the 
two chain regime again. But this time there is no An- 
derson localization at all, and £ = 00. This is the regime 
of Fig. 4. These successive transitions for fixed t± are 
conveniently described as a function of K p on figure 5. 
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Fig. 5. The behavior of the localization length £i oc . as a function 
of K p for fixed t± and D. For K p small enough, the gap in 
fin is well formed and we have the pinned CDW n or two chain 
Anderson Insulator. 9 ' For K p large enough, the gap in 0u is also 
robust to disorder and we obtain the two-chain metal. 9 ' In the 
intermediate regime, disorder can compete with the formation 
of a gap in 8u and establish a single chain regime. 



becomes classical, but now its expectation value in the 
ground state of the system is -y=, making disorder ir- 
relevant. In the pure system, this corresponds to the 
s-wave superconductor 9 ) (SC S ) obtained for attractive 
interactions and there should be no Anderson localiza- 
tion. This is indeed what is obtained since disorder is 
scaled to zero(cf. Fig. 4). We define K" [t - (t± , V) as the 
value of K p above which this regime is obtained. 

3.2 Transition between the different regimes 

We discuss qualitatively the crossovers between the 

different phases in the system. 

We have to compare three correlation lengths (or the 

associated energy scales), the single chain localization 



The crossover from one chain to two chains can be ev- 
idenced easily in the numerical solution of the RG equa- 
tions. In Fig. 6, l/£ioc. as a function of K p is shown 
with fixed D = 1(T 5 and t± = 10 _1 , • • • , lO -5 where 
only the regime for K p > I is shown. For K p — 1 small 
enough, all curves display the same behavior, which is 
the one obtained for t± = 0. In this range of K p , disor- 
der causes Anderson localization at a lengthscale much 
shorter than the correlation length of 6\\ and the system 
is in the single chain regime. When K p — 1 is larger, the 
localization length is increased by the combination of t± 
and the stronger attractive interaction. The change of 
behavior of £i oc . from single chain to ladder occurs ear- 
lier when t± is increased. For K p > K p rlt -(t±), one has 
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£i oc . = oo. We remark that (as could be expected) K" 
is a decreasing function of t± . 
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Fig. 6. The i4" p -dependence of l/£i oc . and 1/£ii (dotted curve) 
with fixed t± = 10 _1 (solid curve) 10 (dotted curve), 
10 -3 (dashed curve), 10~ 4 (dash-dotted curve) and 10 -5 (solid 
curve) and fixed D = 10 . 



It has been well known that in a spinlcss fermion sys- 
tem the energy scale £q oc . = -f 1 - for Anderson local- 
ization is enhanced by repulsive interaction and is sup- 
pressed by attractive interaction. The latter case in the 
two-coupled spinless fermion chains system results in a 
strong competition between Ay and Ei oc .. Actually, by 
comparing An in the absence of the impurity with E\ oc . 
in the absence of the interchain hopping, we obtain the 
boundary between the localization and the derealization 
as (see Appendix B for a derivation) 

(3-2K) 



T> = Co 4t ± exp 



2K -2 



1 



(3.1) 



where Co is a constant. K p :Tlt (t±,T>) can be obtained by 
solving Eq. (3.1) for fixed t±,T>. 

If we compare the transition in the two chain system 
with the one in the single chain system, we note in the 
latter case that the transition is controlled by a single 
critical point at K p = 3/2,1? = 0. As a result, univer- 
sal exponents are obtained at the metal insulator tran- 
sition. This is not anymore the case in the two chain 
system, and we expect non-universal exponents at the 
metal-insulator transition. These exponents might not 
be related in a simple way to if™*- since K crlt - only con- 
trol the exponents associated with charge fluctuations. 
It might be interesting to consider the dependence of 
such exponents of the various correlations with t± and 
V. This will be left for a future study. 

Following a similar line of thought to the one that 
lead to Eq. (3.1), the boundary between the single chain 
localization regime and the two chain localization regime 



is obtained as: 

V = CA 4t ± exp 
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where C\ is another constant. The corresponding phase 
diagram is plotted on figure 7. 
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Fig. 7. The phase diagram predicted from scaling arguments Eq. 
(3.1) and Eq. (3.2). The line separating the two localized regime 
is a crossover line. 



We can compare the predictions of these scaling argu- 
ments with the results of the numerical solution of the 
RG equations (2.8). In Fig. 8, such a comparison is 
shown. The localized region increases with decreasing 
t±. These results are compared with those obtained by 
eq.(3.1) where Co = 0.0207 is chosen so as to fit the data 
for t± = 0.1. It turns out that this formula fits well 
especially for D < t±/10. 
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A plot of V CIlt - vs. K p with the fixed t± 



(solid 
" 4 (dash- 



10" 
curve), 10 — 2 (dotted curve), 10~ 3 (dashed curve), 10~ 
dotted curve), 10~ 5 (solid curve) and 10~ 6 (dotted curve). The 
dots for t± = W^ 1 and 10~ 6 arc calculated from eq.(3.1). 



The effect of t± on K c p rit - 



■ (T>) defined as the value of K p 
below which the system is localized for a given disorder 
strength T> can also be analyzed using the RG equations. 



In Fig. 9, K" as a function of t± is shown for var- 
ious disorder strengths T> = 10~ 7 ...10~ 2 . For t± —* 
0, K c J lt - goes to the single chain value 3/2. As t± is 
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vith the fixed D = 10~ J (solid 
curve), 10 — 3 (dotted curve), 10~ 4 (dashed curve), 10" 5 (dash- 
dotted curve), 10 — 6 (solid curve) and 10~ 7 (dashed curve). 



increased,-ftTp rlt, (£jJ decreases to K" lt * > 1. The pre- 
cise value of this quantity could in principle be obtained 
by setting £j_(0) = oo in Eqs. (2.8). With increasing 
D, K C p lt - increases and the regime for localization is en- 
larged. 



ing from one chain to two chains. 14 > The solid curve 
for K p = 1.1 (attractive interactions) shows a transition 
from an Anderson localized to a superconducting state at 
V = X> crit - ~ 0.12 x 1(T 6 i. e. between Fig. 3 (D = 1CT 6 ) 
and Fig. 4 (T> = 10~ 7 ). Comparing with the single chain 
case t± = 0,K p = 1.1 (dashed curve), we see that well 
above D ~ £> crlt - l/£i oc . is enhanced slightly. Thus, in 
this regime the interchain hopping enhances localization 
even for attractive interactions. This is a consequence 
of a change of sign of yj_ and yf that occurs at large 
lcngthscales , which implies that at large enough length- 
scales the system behaves as if it were repulsive. This 
effect could not be captured by the approximations of 
Ref. 9 since they only apply to infinitesimal disorder. 
Finally, we note that for stronger attractive interactions, 
K p = 1.2, the disorder strength needed to cause a tran- 
sition from the metal to the Anderson insulator is in- 
creased as could be expected. Note however that this 
increase is by 3 orders of magnitude although interac- 
tion strength is only doubled. Thus, we see that the RG 
allows us to recover the qualitative behavior obtained for 
the weak disorder case, 9 ) but also allows us to analyze 
in details the interplay of disorder and interaction. 

§4. Conclusion 



3.3 



We have presented the results of a renormalization 
group approach to the two chain spinless fermion sys- 
tem with disorder. This approach confirms the results of 
behavior of the localization length with disorder strengthRci. 9 in the limit of disorder small enough not to mod- 
ify the gaps of the pure system, namely the strong re- 
inforcement of Anderson localization by interchain hop- 
ping in the case of repulsive interactions and the der- 
ealization by interchain hopping in the case of attractive 
interactions. It also enabled us to obtain new results 
in the regime where disorder and interactions are of the 
same order of magnitude. In the latter regime, a strong 
enough disorder can suppress completely interchain hop- 
ping and maintain the system in a single chain regime. 
In the repulsive regime, disorder, interchain hopping and 
interaction reinforce each other leading in particular to 
an enhancement of the gap to interchain charge excita- 
tions with respect to the pure system. In the attrac- 
tive regime, disorder and interchain hopping compete 
with each other. This leads in the delocalized phase 
to a reduction of the gap to interchain excitations by 
disorder. Finally, an exciting effect is the existence of 
derealization by attractive interactions at small length- 
scales followed by Anderson localization at large lcngth- 
scales. This might lead to a non-monotonic behavior of 
the charge stiffness of a disordered spinless fermion lad- 
der as a function of its length. It might prove worthwhile 
to investigate such an effect numerically Another inter- 
esting subject for future studies is the behavior of critical 
exponents at the localizcd-delocalized transition. 




Fig. 10. The quantity l/£i oc . as a function of T> with the fixed 
K p = 0.8, 0.9, 1.0, 1.1 and 1.2 where t± = 0.1. For the compar- 
ison, l/£ in the limit of small t± is shown by the dashed curve 
for K„ = 0.9 and 1.1. 



In Fig. 10, the quantity l/£] 



loc 



of V is for fixed values of K n 



is plotted as a function 
For repulsive interac- 



tions {Kp < 1), we see that l/£i oc . is always enhanced 
by t± in qualitative agreement with Ref. 9. For K p = 1 
corresponding to the absence of intrachain interaction, 
the effect of t± on l/£i oc . is negligible. This is in agree- 
ment with the fact that for non-interacting fcrmions the 
localization length is only enhanced by a factor 2 by go- 
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Appendix A: Derivation of RG equations 

To derive RG equations, we use an Operator Product 
Expansion (OPE) approach. 15 ' in the absence of dis- 
order, the relevant operator product expansions are the 
following: 

cosv / 80(x,t)cosv / 8</)(x',t')- | [l - 4 [(x - x') 2 {d x <j)) 2 

+ (r-r'fid^) 2 ]] 

i V86(x, t) cos V86(x', t 1 ) ~ - [1 - 4 [(a: - x') 2 {d x 9) 2 

+ (T-T'f(d T ef]\ (A-i) 

To be able to make use of these OPEs, we must express 
d x . T 6 as a function of d XtT (j>. Using the equations on 
motion in imaginary time, we find: 
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/\2 r 



{d T r P ? + (d T 0< 



■cos 



iVS(f) a (x,T) 
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- cos 
2 



V2 [<p a p (x, t) - <j) b p (x, t')] cos V2</^(x, t) cos \/20Jj (x, r') 

xcos\/8^|j(a;",r")~ 
V2 [4>°(x, t) - <j> b p {x, t')} cos V2^(x, t) cos s/2<j>\ (x, r') 

(A-6) 
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cos V2 [<j> a p {x, t) - <j) b p (x, t')] cos V2d a (x, r) cos V2$ b (x, r') 

xcosV86b{x",T") ~ 
- cos V2 [<j) a p {x, t) - ^(a;, r')] cos V2W(x, r) cos V29 b ,(x, r') 
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This leads to the RGE (2.8). 
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This leads to the following final form of the second OPE: 
cos VS9(x, t) cos V8d{x', t) - 

^2 , „.2/_ _/\2/a /)\2 
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Appendix B: Derivation of localization-delocalization 
boundary for the attractive inter- 
action 

In the absence of disorder, the antisymmetric gap is 
given by the coupled RG equations: 



(A3) 



dK \\ 



Using Eq. (5:13) of Ref. 15 we recover the RG equations 
of Refs. 10, 12. 

In the presence of disorder, we have to use replicas 
in order to average over disorder. The replicated action 
(that we use only for the purpose of perturbative RG) is 
the following: 



i(y 2 -A 2 yiJ (4A||t ± )) 
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di " v #ii 
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(na) 2 



To analyze these equations, we remark that if t± — 0, 
y± = yf and K\\ = 1 there is no flow in the RG equations. 
Let us consider the limit £_i_(0) « 1. In such case, we 
can make the approximation: Jo(AK\\t±) ~ 1 for small 
x cosv^^ji {x,t) cosv / 2</'||(a;,r') I. As a result, there is no flow of K\\,yf,y±. Neverthe- 
less, i± flows as t±(l) = i±(0)c l (we have neglected y\ 
compared to 1). The approximation Jo(4Aii£_l) ~ 1 will 
therefore break down at the scale I* when 4A||£j_(0)e' ~ 
1. Beyond the scale I*, we make the assumption that 
Jo(4Knt±) ~ 0. Then, y± disappears from the RG equa- 
tion of An and the only remaining RG equations are: 
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The extra OPEs needed to derive the Renormalization 
Group equations in the presence of disorder are: 
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The initial conditions for these RG equations are with 
;V2[r p (x,r) - </>>,r')] cos V20f (s,r) cos y/20\{x,T>) ~ our approximations: 

A||(D = A||(0) = 1 , 
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Vf(i*) = y/(P) ■ 



(B-3) 



It is convenient to introduce the variable yu in (B-2) 



defined by: K\\ 
be rewritten: 



so that the RG equations can 
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= V\\Vf 



(B-4) 



With initial conditions: y\\(l*) = 0. A convenient pa- 
rameterization of the RG equations (B-4) is : 

2/||(0 = 2//(0)sinh#(/) 

y f (l) = y/(0) cosh $(l) (B-5) 

with initial condition 0(J*) = 0. The system (B-4) is the 
reduced to a single differential equation: 

M 
which leads to: 



y/(0)cosh( 



(B-6) 



2(arctan(c e W) - j) = y f (0)(l - I*) 



(B-7) 



The RG flow is cut when y f (l) ~ 1, i.e. for y f (0)e 9 ^/2 ~ 
1. This leads to I — I* — „ ", n< — 1. The correlation length 

£m is then a 



2!//(0) 

ae' - ' e' , i.e. 



^11 = a exp I 7-t- I 

The resulting gap is then Am = ^j- thus 



Ul//(0)/ 4ef ± (0) 
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An = 4etx cxp 
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As long as we can neglect disorder we will obtain the 
expression of An of the preceding section, in order to 
be able to neglect disorder, we need £i c h. S> £ii . This 
implies the criterion: 



DV(3-2*r) <4e j ± _ * 

2j//(0) 



(B-10) 



Thus, we expect the critical disorder D c to be given 



by: 



^■^(■1|) 



(B-ll) 



where Co is a numerical prefactor. 
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